Non-Abelian bosonization of the frustrated antiferromagnetic 

spin- 1/2 chain 



Dave Allen and David Senechal 
Departement de Physique and Centre de Recherche en Physique du Solide, 
Universite de Sherbrooke, Sherbrooke, Quebec, Canada J1K 2R1. 
CRPS-96-11 [ |cond-mat/9606007p 



(May 1996) 

Abstract 

We study the spin-1/2 chain with nearest neighbor (m) and next-nearest 
neighbor (z^) interactions in the regime K2 3> «i, which is equivalent to two 
chains with a 'zig-zag' interaction. In the continuum limit, this system is 
described in term of two coupled level-1 WZW field theories. We illustrate 
its equivalence with four off-critical Ising models (Majorana fermions). This 
description is used to investigate the opening of a gap as a function of k\ and 
the associated spontaneous breakdown of parity. We calculate the dynamic 
spin structure factor near the wavevectors accessible to the continuum limit. 
We comment on the nonzero string order parameter and show the presence 
of a hidden Z2 x Z2 symmetry via a nonlocal transformation on the micro- 
scopic Hamiltonian. For a ferromagnetic interchain coupling, the model is 
conjectured to be critical, with different velocities for the spin singlet and 
spin triplet excitations. 
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I. INTRODUCTION 



One-dimensional quantum antiferromagnets have peculiar properties (exotic ground 
states, gapped excitations, etc.) which are not accessible to traditional methods like spin- 
wave theory or perturbation theory, but require the use of variational, numerical, or field- 
theoretical approaches. In particular, field-theoretical methods have been used successfully 
to predict the existence of an excitation gap iuJhe spin-1 Heisenberg chain0 and the scal- 
ing behavior of the spin-1/2 Heisenberg chain.™. In the latter case, Witten's non-Abelian 
bosonizationa was used to express the spin-| Heisenberg chain as a Wess-Zumino-Witten 
model perturbed by irrelevant interactions. 

In this work, we apply non-Abelian bosonization to the spin-^ Heisenberg chain with 
nearest-neighbor (NN) coupling K\ and next-nearest-neighbor (NNN) coupling k 2 i n the 
regime k 2 This system may also be viewed as two spin-^ chains coupled with a 

'zig-zag' interaction K\ (see figure). This latter representation makes sense physically, since 
such arrangements of atoms occur frequently (see, for instance, Ref. §). The Hamiltonian 
of this system is 

H = «i • Sj+i + k 2 • S i+2 (1) 

i i 

wherein the spins are indexed consecutively along the zig-zag. This model has recently been 
studied by White and Affleck.B We shall extend the somewhat brief theoretical analysis of 
Ref. ^ and describe the system in terms of four massive real fermions (or Ising models). This 
will allow for an easy calculation of the dynamic spin structure factor S(q,u). We will also 
discuss the occurrence of a string (topological) order parameter and the associated Z 2 x Z 2 
symmetry. We shall take some time to describe in sufficient detail the correspondence 
between the system of two spin-| chains and four real fermions, since this is unfamiliar 
ground for many. 

Let us summarize here the main results of this paper and explain its organization. In 
Sect. H we set up the description of two coupled spin-| chains in terms of four Majorana 
fermions. This implies a quick review of the non-Abelian bosonization of a single spin 
chain (its representation as a level- 1 WZW model). The main result of this section is the 
representation (§,0) of the fields of the level-1 WZW model in terms of four Majorana 
fermions and their associated order and disorder fields. This representation allows for a 
representation of the spin operator and interchain interaction with the help of Eq. @. In 
Sect. [II, we write the interchain interaction (and the marginal intrachain perturbation) in 



terms of the four fermions and discuss the renormalization of the couplings and velocities. A 
mass scale m ~ k 2 exp — (/t 2 //ti) appears dynamically and provides a mass for the fermions, 
accompanied by a spontaneous breaking of parity. Lorentz invariance is explicitly broken by 
the interchain interaction and one of the four fermions acquires a distinct mass and velocity. 
In Sect. [IV| we set up the calculation of the spin structure factor S(q, u) (the imaginary part 
of the dynamic spin susceptibility) near the four wavevectors available to the continuum 
limit: q = 0, n, ±7r/2 (when considering wavevectors, we regard the system ([!]) as a single, 
frustrated chain and not as two coupled chains). The single-spin excitations appear at a 
frequency u ~ m near q = ±7r/2 while a two-particle continuum appears near q = and 
q = ir, like for the spin-ladder. In Sect. |V| we show that the nonlocal string order parameter 
of Ref. [7| is nonzero in the ground state, and how this breaks down a hidden Z 2 x Z 2 symmetry 
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of the system. We also perform an exact (i.e., discrete) nonlocal transformation of (jl]) that 
reveals this Z2 x Z2 symmetry. In Sect. [VJ we discuss the difference between this system and 
the usual spin ladder and address the case of ferromagnetic interchain coupling. A quick, 
largely notational review of WZW models and of the Ising model (Majorana fermion) is 
given in Appendices A and B. 



II. CONTINUUM DESCRIPTION OF TWO SPIN CHAINS 



A. Non-Abelian bosonization 

From the Bethe Ansatz solution we know that the spin-| Heisenberg chain Hamiltonian 

H = K J2Si-Si+i (2) 

i 

is critical. It was also argued by Affleck! that this critical point is well described by a 
level-1 Wess-Zumino-Witten conformal field theorycl (cf. Appendix A). This equivalence is 
demonstrated by starting from the half-filled Hubbard model with hopping integral t and 
on-site repulsion U and taking the continuum limit. The charge degrees of freedom are then 
described by a Bose field tp which becomes massive for arbitrary small U, while the spin 
degrees of freedom are described by the level-1 WZW model. At U = the characteristic 
velocity v of the WZW model is simply the Fermi velocity Vp = |t|oo ( a o is the lattice 
spacing). For U > the velocity v of the spin degrees of freedom is renormalized by U and 
differs from the velocity of the charge excitations (spin-charge separation). Moreover, the 
continuum limit of the Hubbard Lagrangian contains an additional term: 

d = -XJ a J a (3) 

where J a and J a (a = 1,2,3) are the left and right components of the SU(2) currents of 
the level-1 WZW model and A ~ U/\t\ (we will work in the imaginary-time Lagrangian 
formalism; C denotes the Lagrangian density). This perturbation is marginally irrelevant. 
Thus, at lone, distances, the spin degrees of freedom are exactly described by the level-1 
WZW model! 

Additional perturbations to the Heisenberg Hamiltonian (Q) may be expressed in terms 
of WZW fields by using the following continuum-limit expression for the spin operator SjU 

— = ^ (J a W + ■/*(*)) + (-l) x/ao eTr(g(x)r a ) (4) 

CLq Z7T 

where r a are the usual Pauli matrices and g(x) the fundamental WZW field (an SU(2) 
matrix). The factor (— Y) x l a ° alternates from one site to the next and G is a nonuniversal 
constant. The first two terms of Eq. (|j) constitute the local magnetization and the last 
term is the local staggered magnetization. 

Let us now turn our attention to the system (|TJ). In the regime K\ k 2 and in the 
continuum limit, it may be regarded as two level-1 WZW models, plus some perturbations. 
Let J a and J a denote the SU{2) currents on one chain and J' a and J' a the corresponding 
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currents on the other chain. The first perturbation is marginally irrelevant and given by two 
copies of (Q): 

C 2 = -2\ 2 (J a J a + J' a J' a ) (5) 

where A2 ~ U/\t\. The second perturbation is the interchain interaction (ki). In the 
continuum limit and using Eq. (Q), it can be shown without difficulty to be 

d = 2Ai(J a + J a ){J' a + J' a ) (6) 

where Ai ~ Ki/\t\ > 0. 

The relevance or irrelevance of a perturbation is determined, as a first approximation, 
from the scaling dimensions of the various fields at the WZW fixed point. In terms the 
conformal dimensions (h, h) appearing in (|A8|), the scaling dimension is A = h + h and 
the planar spin is h — h. Since the conformal dimensions of J a and J a are respectively 
(1, 0) and (0, 1), a perturbation of the form J a J a (like ([5])) is marginal, while a perturbation 
of the form J a J a violates Lorentz (or rotation) invariance. In fact, it renormalizes the 
characteristic velocity of the theory (see below). The interaction (||) is marginal, except for 
a velocity renormalization. 

WZW models, although they possess conformal invariance, are not easy to deal with, 
especially in what regards the calculation of correlation functions. In some cases (i.e., for 
some values of the level k) the WZW model is equivalent to a theory of free fields. Then the 
calculation of correlation functions becomes an almost trivial task and the overall analysis 
is much simplified, in particular the study of the vicinity of the fixed point. Such a free-field 
description is possible in the case of two coupled level- 1 WZW models: two such models are 
equivalent to one level-2 WZW model, plus one Ising model (or real fermion, see appendix 
B). This equivalence was already used in Ref. || to describe the spin ladder with .bond 
alternation. Moreover, the level-2 WZW model is equivalent to three Ising modelsB We 
thus have three different ways of describing the system (TJ) in the continuum limit: 

WZW fc= i ® WZW fc= i (7a) 
WZW fc=2 <g> Ising (7b) 
(Ising) 4 (7c) 

The representation (0b) may be useful from the point of view of symmetry since the interact- 
ing terms break down the SU(2) x SU(2) symmetry to SU(2). However, the representation 
(0c) is more practical for calculations since it is made entirely of free fields and its off-critical 
(/«i 7^ 0) behavior may be characterized by ordinary fermion mass terms. All is not trivial, 
however, since the Ising model contains order and disorder fields in addition to a real fermion 
field, and these three fields cannot be expressed locally in terms of each other. Nevertheless, 
their correlation functions are known. An additional difficulty comes from the breaking of 
Lorentz invariance by the perturbation ([]). 

We identify operators in two different languages of (0) by requiring their operator product 
expansions (OPE) to be compatible. The OPE for the WZW models and the Ising model 
are given in appendices A and B. The correspondence of operators belonging to the pictures 
(0a) and (0c) is the object of the next subsection. 
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B. Description in terms of four Ising models 



The WZWfc=i model cannot be simply represented in terms of two Majorana fermions, 
even if the central charge is the same in both cases (c = 1). The reason is the nonexistence 
of a real, spin-^ representation of SU(2). However, two copies of WZWfc=i is equivalent 
to an 5*0(4) WZW model, and the latter group admits a representation in terms of four 
real fermions. A representation of the WZW currents J and J' in terms of four Majorana 
fermions ^1,2,3,0 follows immediately and its structure bears a strong resemblance with the 
chiral generators of the Lorentz group: 



J x = 




- ^2^3 


J 2 = 


|i(-0 2 -0o 


- IpZ^l 


J 3 = 







(8a) 
(8b) 
(8c) 



(the corresponding expressions for J' a are obtained by reversing the sign of i/jq). Using the 
OPE's ( |B1| ) and Wick's theorem, it is a simple task to check that the OPE's ( |A12| ) are 
satisfied. 

A representation of the matrix fields g and g' (the staggered magnetizations of the two 
chains) in terms of Ising fields is also needed in order to calculate correlation functions, and 
may be found in the following fashion. First, since g and g' have conformal dimensions 
(i, |), they must be products of four order and disorder fields, such as 01O2O3O0, cri/i 2 03/io, 
and so on (there are 2 4 = 16 such products). Second, the action of each of the currents J a , 
J a , J' a , J' a may be calculated on these 16 products, using the OPE's (|B3|e-h). The result is 
a 16- dimensional matrix representation of the currents. According to the OPE ( |A11| ), the 

field 011 is an eigenvector of J 3 with eigenvalue — ~. Once such an eigenvector is found, one 

22' ' 2 _ 

may apply on it the other components of the currents J and J and thus obtain the other 
components of g. Only one eigenvector allows a nontrivial solution (i.e., nonzero values of 
all the components of g). The same procedure is used for g', with the currents J', J'. At 
last, one finds the following representation (we used the decomposition (|A15|) ): 



go = OW2OW0 - ^1^2/^3/4) 

91 = /UlC2C3A*0 - O"lA i 2/i3O"0 

92 = O"D"2^3A*0 + AHO"2/U 3 er 

93 = 0102/^0 - ^l/^OWo 



g' 



92 



01 0-2(73 a + /il/!2/i 3 /i0 



-O-lyU2O"3yU0 + ^2^0 



-OiO 2/ U 3/ U - yUi/i2030o 



(9) 



Note that the OPE's J a (z)g' i (w,iD) ~ and J' a (z)gi(w,w) ~ are satisfied, as they should: 
the two chains are independent at this stage. 

It is also possible to calculate the OPE of g with itself, with the help of Eqs (|E3|a-d) . 
This is a bit tricky, since one must remember to anticommute the different disorder fields. 
With the normalization chosen above and omitting terms that do not diverge as z — > w, the 
end result coincides with Eq. ( |A16|) for g and g', plus the OPE g^z, z)g'j(w, w) ~ 0. Thus, 
the representation (||,||) is a complete and faithful representation of two independent copies 
of the WZW fc= i model. 



5 



III. VELOCITY RENORMALIZATION AND RG ANALYSIS 



We are now able to write down the Lagrangian associated to the continuum limit of ([!]) 
solely in terms of real fermions. The noninteracting part Co, equivalent to the two level-1 
WZW models, is the free-fermion Lagrangian: 



1 3 

i=0 



(10) 



where Vq = • • ■ = v 3 = v is the velocity of spin excitation in isolated chains. The 2tt factor 
in fllOD is needed for consistency with the OPE (pT|). 

The interacting terms (fUD may be expressed in terms of the following operators: 



2 = ^oV^l^l + ^2^2 + ^3^3) 



The interaction (ffl) is simply 



A = -A 2 (0 1 + 2 ) 



(11a) 
(lib) 



(12) 



The translation of (|(|) requires more care, however, since it implies regularized products of 
identical fermions. The following OPE must be used to extract the regular terms: 



z — w 

1 

z — w 



+ 2{z-w)T^\w) 



+ 2(z -w)T {i \w) + 



(13) 



where (T«,T«) is the energy-momentum tensor ( |B2|) of the i-th Ising model. The result is 



C 1 = \ 1 (0 1 -0 2 )+\ 1 



i=l 



(14) 



Apart from the energy-momentum terms, this interaction coincides with the marginal in- 
terchain interaction obtained by Shelton and al.0 using Abelian bosonization. The effect of 
the energy-momentum tensor is simply to renormalize the speeds V{ of the fermions. Indeed, 
consider the Lagrangian (we have restored the velocity v in the interaction term) 



v - - - 1 

£ = — (ipdip + ipdtp) v\(ipdip + ipdif)) 

2n 2 



(15) 



(cf. Eq. ( P2| )) where A is a dimensionless parameter. One may combine the energy- 
momentum tensor with the kinetic term and this amounts to the following renormalizations 
of the speed and fields: 



ttA 



V — > V- 



1 



1 + ttA vr ' r/ VI + vrA 

In the case at hand, the velocity vq is renormalized differently from -i^, t> 2 , t^: 



(16) 
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1 + 371 " Al 1 -7rAi 

u o ^o^ 5— r- «i -> Ui^— r- = 1, 2, 3) (17) 

After the field renormalization, the interaction Lagrangian takes the following form: 

£int - = (l + nxly 01 - (l + nX^-SnX,) 02 (18) 



The 0(4) symmetry of the fixed-point, obvious in the Lagrangian (10), is violated by the 
interchain coupling Ai, both in the interaction (|T8| ) and by the distinct renormalization of 
v . 

The interaction terms ( |IS| ) are marginal, since they have conformal dimensions (1,1). 
Their behavior under renormalizat ion-group flow is characterized by their /3-functions. In- 
stead of calculating the latter in the usual way (a one- loop Feynman diagram), let us follow 
Polyakov,EJ who has shown that the /3-functions of a critical system perturbed by marginal 
terms are related to the coefficients of the operator algebra. Explicitly, consider the per- 
turbed action 

S = So + J^Xi f d 2 x <j>i(x) (19) 

i 

where S is the fixed-point action and the <f>i(x) are marginal operators (h = h = 1). Let 
the operator algebra be of the form: 

mm - ^# (20) 

\x y\ 

Then the renormalization-group flow of the couplings Ai, characterized by the /3-functions 
flijk(L), is (in Euclidian space-time) 

Pijk^j^k = —ZTfCijkXjXk (21) 



dlnL 

If we apply this method for the perturbation (|i"8|), we must use the OPE 

~ -rj^Tp (22) 

and realize that the eigenmodes of the RG flow are the operators 

JC+ = O x + 2 /C_ = Oi - 2 (23) 

which have the OPE 

JC ± (z)lC ± (w) ~ — ^—rj - -J—)C ± (w) + —L—0(T + f ) + ■ ■ • 

\z — w\ 4 \z — w\ 2 \z — w\ 2 

K + {z)K„{w) ~ . 1 n(T + f) + --- (24) 
\z — w\ z 

where 0(T + T) stands for terms containing the energy-momentum tensor which, although 
they have the right scaling dimensions, also have nonzero planar spin and do not contribute 
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to the beta functions. The terms in \z — w\~ A only contribute to a shift of the vacuum energy 
and will be ignored. The interaction ( |T8D may be expressed as a linear combination of the 
operators /C±: 



with 



Ant. = A + /C + + A_/C_ (25) 



\ 1 fAi + A 2 , Ai — A 2 , 

x± ~ 2(1 + ttAO 1 t T^a7 + TT^ t (26) 



The RG equations obtained from (|2lD and (|J) are 

c/A+ 



dlnL 



ttA^ (27) 



If Ki <C «2 the starting point is Ai 1 and A2 ~ 1, thus A + is negative and A_ positive 
and small. Under this flow, A + renormalizes to zero (it is marginally irrelevant, like for an 
isolated spin chain) and A_ is marginally relevant. By following the RG flow until L = £ (the 
correlation length), we conclude that £ ~ exp(l/A_): a dynamical length scale £ has set in. 
If we concentrate on the Heisenberg model ([TJ) without referring to the underlying Hubbard 
model, the characteristic spin velocity should be v ~ k^o and the dimensionless coupling 
constant A_ should be ki/k 2 . Thus the dynamical length scale is £ ~ exp— k 2 /«i- This 
conclusion was reached in Ref. ^in the (WZW^i) 2 representation. We shall conclude from 
this that the fermions have acquired a mass ~ v£~ l . The first three fermions (^1,2,3) 
have a common mass m, while ipo has a slightly higher mass itiq (in absolute value), since 
Vq > v\ = v% = V3. An additional velocity renormalization will take place during the RG 
flow, but cannot be calculated by the above method. In a diagrammatic technique, it would 
show up at two loops, in a self-energy correction. The overall velocity renormalization is 
important, since without it the long-distance theory would have an 0(4) symmetry and all 
four fermions would have the same mass (up to a sign). 

In order to clarify the significance of these dynamically generated mass scales, let us 
consider the following model: 

I N 1 

£ = — v (tpidtpi + ipidipi) + -A ipiipiijjijj (28) 

Z7T z — ' Z * — ' 

1=1 tyj 



This model would be equivalent to the Lagrangian (|TU|,|TB]) if all velocities were equal, if 



A + = and if N — 4, with i/j 4 = ipo, 4>4 — ~4'o (Kramers- Wannier duality). The model 
has 0(N) symmetry and a mass gap arises non-perturbatively in the spectrum if A > 0. 
To see this in a mean-field approach, we assume that (ipiipi) = is 7^ (no sum over i) and 
determine e self-consistently. Let us make the substitution 

^i^i — ► ie + (29) 

in the Lagrangian, neglecting terms quartic in if), which is equivalent to a large- N approxi- 
mation. We find the Lagrangian of N massive fermions: 

N 

C = H{^ + Wi) + *A(iV - l)eVi&} (30) 
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where the mass is m = 27r\(N — l)e. This mass may be determined self-consistently, using 
the following expression for the Green's function of real fermions:@ 



/Mh, -ikx 
7T K z + m z 

— (31b) 

M0Mx)> = -i\m J *± (31c) 

where x and k stand respectively for (vt,x) and (iuj/v, k). The mass m is determined from 
the self-consistency condition 

m f d 2 k e~ lkx . . 

m / — TJT-ZT, (32) 



27rA(JV-l) J n k 2 + m 2 
whose solution is, besides m — 0, 

- = ±.Aexp-^l__ (33) 
where A is a momentum cutoff. This solution exists only for positive A. 



Returning to the Lagrangian ( |10| , |18l) with A + renormalized to zero, all velocities equal 
and A = A_, this implies a mass gap m ~ v Aexp(— l/67rA_), or m ~ /t2exp(— k^/Qitki) if 
the characteristic velocity (of order k 2 ) is restored. Since ('04, ^4) = (4>o,— 0o), the mass 
m of the fourth Ising model is equal to — m, if v = V; L . Since v > Vi, we conclude that 
— m > m. 

A short remark about the sign of the mass: from the Ising model viewpoint, this sign 
simply indicates on which side of the transition we stand: By convention, m > in the 
disordered phase ((/x) 7^ 0) and m < in the ordered phase ((a) 7^ 0). Of course, it is the 
absolute value \m\ which occurs in the dispersion relation of the fermions. 

The appearance of fermion mass terms breaks the diagonal Z2 symmetry (-0j,0j) — ► 
(tpi, —ipi) (i = — 4) of the full Lagrangian ( |l0| , |nj|) . Thus, the ground state must be doubly 



degenerate and the condensate (ipiipi) 7^ picks one of these ground states, the theory of 
massive fermions describing excitations above that ground state only. This is consistent 
with the Lieb-Schultz-Mattis theorem, which states that a half-integer spin chain with local 
interactions and no explicit parity breaking has either no gap, or else has degenerate ground 
states. 



IV. SPIN STRUCTURE FACTOR 

In a recent paper, Rao and SenEl have argued that dimerized spin chains with second 
nearest-neighbor interactions admit possibly three different phases (here we use the word 
'phase' to distinguish regions where the spin structure factor S(q) is not peaked at the same 
value of q). They name the three phases as follows:i a Neel phase (S(q) is peaked at tt), 
a spiral phase {S(q) is peaked at an intermediate momentum between 7r and 7r/2) and a 
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collinear phase (S(q) is peaked at 7r/2). In view of the numerical results from Chitra and 
al.@, the collinear phase should not be stable for the spin-| chain. At first sight, there are 
two paths that S(q) may follow to go from the Neel phase to the collinear phase. The first 
possibility is for the peak of S(q) to move continuously from n to tt/2, thus going through 
a spiral phase. The second possibility is for the peak of S(q) at 7r to progressively decrease 
in amplitude while a second peak at ir/2 progressively appears; the system might then go 
through a dimerized state. In view of this, the question of the existence or not of a spiral 
phase for the frustrated spin-^ chain arises. To answer this question, we need to known 
how the spin structure factor evolves as a function of the ratio K\j k^. Unfortunately, in the 
present continuum approach we can only calculate the spin-spin correlation function near 
q = 0, q = tt and q = ±7r/2. As seen below, this calculation is also interesting from the point 
of view of symmetry and allows to relate the elementary spin excitations to the fermions fy. 

The main conclusion of Sect. |l|is that the system (H) may be described in the continuum 
limit by four noninteracting real fermions: three with mass m > and velocity v, and one 
with mass < —m and velocity v > v. The spin operator S« is represented in terms of 
these fermions through the relations (|4]), (§) and (|9]). The z-component of the spin density 
has the following expression near the wavevectors accessible to the continuum limit: 



Sg^ oc 4>i?p2 + i>i*h 

OC 1p 1p 3 + ^0^3 
Sq~v/2 °C ^2^0 
Sq tt/2 0(1 Vl^l^HO 

Thus, the spin-spin correlation function near q = takes the form 
X i0 \x, T) (X ( 0iV>2 + V^ 2 ) (x, t) 0iV>2 + M2) (0, 0)) 
while near q = it, it takes the following form: 

X (T \x, t) cc ( (V^o + ^0) (x, t) (V^o + ^0) (0, 0)) 



(34) 



(35) 



(36) 



In the first case (q near 0), the two fermions have the same mass and velocity, while in 
the second case (q near 7r) they have different masses and velocities. Consider the case 
of two fermions with different masses (m and m!) but identical velocities (for simplicity). 
The imaginary part of the Fourier transform of the spin-spin correlation function - i.e., the 
imaginary part of the dynamic susceptibility, or the spin structure factor S(q, uj) 
calculated from the propagators fl3"T|): 



u 



uj- 



S(q, u>) oc — —(m + m) H — -(m — m) 



,x2 u 2 + q 2 



>\2 



(m + m) (m — m!) 



may be 
(37) 



where u and s are defined by 



f 2 , 2 I2\2 a 2 2 

[s +m — m ) — 4m s 



2 22 
uj — v q 



(38) 
(39) 



(in this expression we have returned to real frequencies). In the case m = m' this result 
becomes the spin structure factor near q = 0: 
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SW(q,u)oc , (40) 
s 3 v s 2 — 4m 2 

Neglecting velocity renormalization (vq = Vi), the model is 0(4) symmetric at long distances 
and m Q = —m. Then the dynamical susceptibility near q = tt would be given by the more 
general expression (|3~7|) with ml = —m. The expression of S^ n \q,uj) appropriate for the 
more realistic case vq ^ Vi can be obtained in closed form but is too cumbersome to display 
here. 

The fact that S , ^°- ) (0, uS) = reflects the conservation of the total magnetization. This is 
not the case for the total magnetization of each chain, since m ^ m', unless the two chains 
are decoupled (m = m' = 0). Thus, SW(Q,u) ^ 0. 

According to Eq. Q53p, the magnetic susceptibility near q = ±tt/2 is a product of four 
two-point functions of the Ising model, involving order and disorder fields. For instance, 
near q = tt/2, 

X {n/2) (x,r) oc ((jiifi2<r3<7o)(0,0)(jJ>it*2<r3<ro)(x,T)) 

= C 2 (mr)C(mr)C(m r) (41) 

where C(R) and C(R) are respectively the two-point functions of the order field and disorder 
field, as a function of the reduced distance R = \/x 2 + v 2 t 2 /^ = mr. These functions are 
knownEj and their leading asymptotic behavior is, in the disordered phase of the Ising model, 

C{R) = -K (R) + 0(e- 3R ) C(R) = a{i + + 0(e' m )\ (42) 

7T [ Q7TK J 

where A is some constant and _Ko,i are the modified Bessel functions. If the argument of C 
is negative (i.e. for a negative mass), we perform a Kramers- Wannier duality transformation 
and identify C(—R) with C(R). Thus, the leading asymptotic behavior of the susceptibility 
near q = tt/2 is 

X [n/2) {x, t) oc C 2 (mr)C(\m \r)C{mr) oc K {mr) + 0(e" 2mr ) (43) 

Notice that K (R) is the real-space propagator of a free boson of mass m. Thus, its Fourier 
transform is ~ (k 2 + m 2 )~ l . Going back to real frequencies, the imaginary part of the 
susceptibility has a pole at u = \J (vk) 2 + m 2 , plus an incoherent part starting at uj = 2m: 

S^l 2 \q,uu) oc - — -5{uo — a/ {vk) 2 + m 2 ) + incoherent part (44) 
\oj\ 

The magnetic susceptibility near —tt/2 is obtained by Kramers- Wannier duality: 

X ( ~" /2) (£,0 oc C(mr)C 2 (mr)C(\m \r) oc K (mr) 2 K (\m \r) + 0(e- 2mr ) (45) 

The associated spin structure factor has no single-particle peak, but instead a continuum 
starting at uj = 2m + \m \. Thus, the single-particle magnetic excitations live around 
k = tt/2, whereas the excitations near k = 0, tt have a two-particle behavior and those near 
k = —tt/2 have a three-particle behavior. 

The above analysis assumed that the mass m was positive, which amounts to choosing one 
of the degenerate ground states, characterized by a short-range order around k = tt/2. If the 
other ground state were chosen, the mass m would be negative and the above analysis could 
be repeated by interchanging the roles of k = tt/2 and k = —tt/2. Thus, the spontaneous 
breakdown of parity is reflected in the nonequivalence of S(k,u) and S(—k,u). However, 
this would be unobservable in practice because of domain effects. 
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V. STRING ORDER PARAMETER AND Z 2 x Z 2 SYMMETRY 



Kohmoto and Tasaki have shownli^l that, for a spin-| chain with dimerization, a string 
order parameter may be defined as den Nijs and Rommelse have previously done for the 
spin-1 chain.i A nonlocal unitary transformation is introduced to show that the nonzero 
value of this string order parameter is related to the breakdown of a Z 2 x Z 2 symmetry. 
More recently, Shelton and al.0 showed that the string order parameter: 

O z (n, m) = exp J m ( S j + Sj) 1 (46) 

I j=n ) 

becomes, in the continuum limit, 

lim (0*(s, y)> ~ K) 2 (a 2 ) 2 - {^) 2 {^f (47) 

\x-y\->oo 

It was also argued in Ref. |H] that the nonzero value of this order parameter is related to the 
breakdown of a Z 2 x Z 2 symmetry which, in the continuum limit, is given by the invariance 
under sign inversion of both chiral components of each Majorana spinor: ip a — > —ip a arid 
ip a — > —ip a (a— 1,2). This must be accompanied by an inversion of both order and disorder 
fields: a a ^ -a a and /i a -> -fi a . 

Here, we expect a nonzero value of the string order parameter for two reasons. First, since 
the SU(2) symmetry cannot be spontaneously broken - according to the Mermin- Wagner 
theorem - the mass of the first Ising model (mi) must be the same as that of the second 
Ising model (m 2 ), in order to keep the symmetry under the exchange of the labels 1,2 and 
3. It naturally implies that they must have the same sign. So, if (ai) ^ (mi > 0) then 
(a 2 ) 7^ (m 2 > 0). Similarly, if \i\ ^ (mi < 0) then (/i 2 ) ^ (m 2 < 0). Secondly, since 
a gap open by the introduction of the interchain coupling, the masses mi and m 2 must be 
nonzero. 

We can also reveal the presence of the Z 2 x Z 2 symmetry without going to the continuum 
limit, i.e., directly from the Hamiltonian (|T|). The unitary transformation U introduced in 
Ref. Ill] consists of many transformations applied in succession. Explicitly, we have: 



U = (D T Y l RDG (48) 
where G performs a rotation of tt about the y-axis on some of the spins: 

L/2 



G = Q£exp 



?7T 

2 l°4j-l ' 



(49) 



I) is a duality transformation (see appendix A of reference^) which introduces intersite 
spins. It is followed by a translation: 

R:r^-(r + -) (50) 



2 V 2, 

The spin on integer sites will be noted a and those on the half-odd integer sites will be 
noted r. The final operation is to make an inverse duality transformation for the r spins. 
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We refer the reader to the work of Kohmoto and Tasak£9 for a full description of this unitary 
transformation. If we apply this transformation to the Hamiltonian (|l|), we find: 

TV 

TTUTT — 1 \ 11 I i _:r_£ _z_z _;g „_z „z^z^z „z 1 

UHU 2^\'' ; + °i - 3 T J - a 3+l - T J 7 J I I - °3 T 3 7 >+l f7 i+l] 

3=1 

+k 2 [<t*t*<t* +1 + a]rjr^ +1 + r/^r;., 

•^•i^-i + ^^/"f-i^-i + } ( 51 ) 

The ex's and r's are sets of Pauli matrices. The new Hamiltonian H = UHU~ l is clearly 
invariant under a rotation of 7r about the x axis applied to the cx-spins alone or the r-spins 
alone. An four-fold degeneracy of the ground state of H in the thermodynamic limit does 
not follow from this broken symmetry since this is not a local symmetry.li3 



VI. DISCUSSION 

The crucial difference between the Hamiltonian (|J) and that of the more familiar spin 
ladder is the occurrence, in the latter, of an interaction term of the form 

Aadder = ^ (gr a )Tr (g'r a ) (52) 

(77 is some constant proportional to the interchain coupling). Since the matrix field g has 
conformal dimensions (7, |), the above perturbation has scaling dimension 1: it is relevant. 
If the representation (Q) and the OPE's (|B3| a-d) are used to express this interaction in terms 
of fermions, one finds 

Tj - - - - 

^ladder = + ^2^2 + ^3^3 ~ 3lpo4>o) (53) 

Z7T 

This coincides with the conclusions of Ref. [10], obtained by Abelian bosonization. The 
mass terms now appear explicitly, with a triplet of mass 77 and a singlet of mass — 377. 
The interchain coupling explicitly breaks the invariance under parity that is spontaneously 
broken in the 'zig-zag' case. If the two rungs of the zig-zag had different couplings (k\ and 
an interaction like ( |52"D would be generated and the gap would have a linear dependence 
on the interchain coupling. Of course, the marginal interaction ( flq) is always present and 
provides an additional renormalization of the masses. As k[ — > k%, the dependence of the 
gap on the interchain coupling should become more and more exponential because of this 
renor malizat ion . 

We were concerned in this work with the regime K\ <C k 2 and the conclusions are 
nominally valid only in this regime, although we expect them to be qualitatively correct 
even for k± ~ k 2 . However, in the opposite regime (acj 3> k 2 ) the system should be treated 
as a single chain and we should perturb around a single WZWfc=i model. This is explained 
in Ref. |3|. The conclusion is that the perturbation is marginally irrelevant if the ratio k 2 /k\ 
is smaller than some critical value, and leads to an exponential gap above that critical 
value. In that regime the ground state is spontaneously dimerized (spontaneous breaking of 
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parity). This conclusion is also valid in the regime K\ <C k^. Indeed, the order parameter 
for dimerization 

d = (S2i ■ S2i-i — S2i ■ S2i+i) (54) 

coincides, in the continuum limit, with the ladder perturbation (|52"D, up to terms that have 
a vanishing expectation value. Translated in terms of the bare interaction couplings X± of 
Eq. ( p5[) and of the masses m, mo, the spontaneous dimerization becomes 

d oc m [m(A_ — A + ) + m (A_ + A + )] (55) 

This is generically nonzero. 

So far we have supposed that K\ is positive, corresponding to an antiferromagnetic in- 
terchain coupling. The ferromagnetic case may be treated just as well. In that case, both 
interaction constants A± of Eq. ( p5|) are negative and thus renormalize to zero: the model 
is equivalent to a theory of four free Majorana fermions, with different velocities. From 
Eq. flT?p with negative Ax, we expect the velocity vq of ipo to be smaller than the velocity v 
°f "01,2,3- Thus, we conjecture that the ferromagnetic model is critical, albeit with two sectors 
having different velocities: a triplet sector equivalent to the WZWfc=2 theory and a singlet 
sector with a smaller velocity. This is not the same as saying that the two chains decouple 
at long distances, since in that case the structure of excitations would be different. This 
conjecture might be tested by exact diagonalizations on finite systems and some information 
on the velocity renormalization might be extracted this way. 
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APPENDIX A: WZW MODELS 

A systematic review of WZW models cannot take place in a regular paper. Here we 
simply recall basic concepts and a few definitions, in order to fix the notation and the 
normalization used in this work. We follow in this respect Ref. II? . 



Wess-Zumino-Witten (WZW) models are defined in terms of a matrix-valued field g 
belonging to a unitary representation of SU(2) (more generally, of a Lie group q) with the 
following action:lS 

S = ^jd 2 x Tr'(d»g-%g) - ^- J tfy e^Ti' {g- l d,gg~ l d u gg- l d p g) (Al) 
where the trace Tr' is proportional to the usual trace operation: 

Tr' = — Tr x s = -s{s + l)(2s + 1) (A2) 

Xa 3 



14 



(s is the spin of the representation). A; is a positive integer called the level of the WZW 
model. The first term of is the usual nonlinear sigma model. The second term is 
topological and is integrated on a three-dimensional manifold B of which two-dimensional 
space-time is the boundary. Its value is independent of the precise form of B (modulo 2n), 
provided k is an integer. 

The fundamental property of the WZW model - enforced by the relative normalization 
of the two terms of the action ( |A1| ) - is its full conformal symmetry. For this reason, it 
is best described in the language of conformal field theory, with holomorphic (or left) and 
antiholomorphic (or right) coordinates 

z = — i(x — vt)= vt — ix 

z = i(x + vt) = vt + ix (A3) 

where r = it is the Euclidian time and v is the characteristic velocity of the model, implicit 
in the covariant notation of Eq. ( |A1| ). The left and right derivatives are commonly used: 

The WZW model has SU(2) symmetry and this entails the existence of a conserved current 
J M , expressed here in its left (z) and right (z) components: 

J = J z = dgg- 1 J=J- z = g- x dg (A5) 

Closely related to its conformal symmetry is the separate conservation of the left and right 
currents: dJ = and dJ = (the SU(2) symmetry is enlarged to a chiral symmetry 
SU '(2) l <8> SU (2) r) . Hence J(z) depends only on z and J(z) on the z. These matrix currents 
may be decomposed along a basis of spin-s generators. For spin-|, we choose 

J( z ) = J a ( z )r a J(z) = J a (z)r a (A6) 

where the r a are the usual Pauli matrices. 

In practice, the action ( |Al| ) is not useful for practical calculations. The traceless, sym- 
metric energy-momentum tensor, which generates local conformal transformations (in par- 
ticular space-time translations) is more useful. Its two nonzero components are given by the 
so-called Sugawara form: 

T(z) = - 1 (J a J a ) f(z) = „ 1 A J a J a ) (A7) 

The notation (...) above stands for a normal ordering (regularized product). The dynamics 
of the theory is determined by the short-distance product (operator-product expansion, or 
OPE) of the various fields. The OPE of T(z) with a local scaling (or primary) field <p(w, w) 
reflects the conformal (or scaling) properties of that field: 

nmw) ~ + d w <p{w,w) 

(z — W) 2 z — w 

where h is the conformal dimension of the field and the symbol ~ means an equality 
modulo terms which are regular as z — > w. A similar expression holds for T and the sum 
h + h is the usual scaling dimension. The OPE of T with itself is slightly different: 
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T{Z)T{W) ~ + - + A9 

The constant c in the most singular term is the central charge of the conformal theory and 
measures the number of degrees of freedom of the theory; its value in the SU(2) WZW 
model is 

c=^ (A10) 

The OPE of the currents J and J with a local matrix field fields reflects its transfor- 
mation properties under the action of SU(2): 

to/ x / \ lr a g(w,w) - n . . . . lg(w,w)r a 
J a (z)g(w } w) ~ -- J a (z)g(w,w) ~ - yy All 

2 2; — w 2 z — w 

The OPE of the current with itself constitutes the so-called current algebra: 

r6 , , (k/2)5 ab , ^ J» 



[z — w) 2 z — w 

lat-\Tbt-\ (k/2)5 ab . J c (w) 

[z — w) z z — w 

J a {z)J b {w) ~ (A12) 

The WZW at level k contains several matrix-valued scaling fields, one for each value of 
the spin s up to (and including) s = k/2. The conformal dimensions h and h of the spin-s 
field are 

h s = h s = (A13) 

The OPE of the various matrix fields is governed by the rule of addition of angular momenta 
and by the constraint that no field of spin s > k/2 occurs in the operator algebra. These 
OPE's were calculated in Ref. |9|. We shall only be concerned with the simplest case (k — 1). 

The level- 1 WZW model has central charge c = 1 and contains a single matrix field g nn 
(n, n — ±2) °f conformal dimensions (|, |). With the proper normalization, its OPE is 

9nn{z,z)g min (w,w) ~ TSnm^nrh (A14) 

\z — w\ 



where e nm is the antisymmetric symbol. We may use the following decomposition: 

3 

2 



1 3 

g(z, z) = -Y J 9a{z, z)r a g a = Tr(r a g) (A15) 



i=0 

The OPE of g with itself is then 



g a (z, z)g b (w, w) ~ ^^(-1)^+ 1 (A16) 
\z — w\ 
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APPENDIX B: THE ISING MODEL 



It is well known that the two-dimensional Ising model is equivalent to a real (Majorana) 
fermion in 1+1 dimension. The critical point of the Ising model corresponds to the massless 
point of the fermion theory and constitutes one of the simplest conformal field theories, of 
central charge c = |. This theory contains a two-component fermion (ip(z),ijj(z)). The 
holomorphic field ip has conformal dimensions (|,0) while its antiholomorphic counterpart 
■0 has conformal dimensions (0, |). Their OPE is 

ib(z)ib(w) ~ 

z — w 

^(,z)^{w) ~ - 

z — w 

tp(z)il>(w) ~ (Bl) 

The product e = iipip has conformal dimensions (|, |) and is called the energy field; it is the 
mass term that takes the model away from its critical point. The energy- momentum tensor 
of the fermion theory is 

T(z) = ~i>d^> f(z) = ~m (B2) 

The critical Ising model also contains an order field a(z, z) which is the continuum limit 
of the Ising spin. This field has conformal dimensions (jq,jq) and is not locally related to the 
fermion field. Indeed, in the transfer-matrix description of the 2D Ising model, the fermion 
field is introduced by a (nonlocal) Wigner- Jordan transformation. The Kramers- Wannier 
duality transformation of the Ising model maps the order field a into a disorder field /x 
which has essentially the same properties, except that (a) ^ 0, (/i) = in the ordered phase 
and (o") = 0, (fi) ^ in the disordered phase. At the (massless) critical point, both fields 
have a vanishing expectation value. All three fields (ip, a, fj,) are mutually nonlocal, which is 
reflected in their OPE by the existence of branch cuts. These OPE's are given below: 

1 1 

a(z, z)a(w, w) ~ - ^ji + 2 \ z ~ w \ 3 ^ £ i w > (B3a) 



\z — w 



/i(z, z)/i(w, w) ~ - — ~\\ z ~ w \ 3/4 z(w, w) (B3b) 



\z — w 



, \ , \ -i{z - w) 1/2 i)(w) +Pi*(z- w) 1/2 ib(w) 
a(z : z)pt(w : w) ~ ii } "\ ~ , /4 (B3c) 

7 *(z- w) l / 2 ip{w) w) l / 2 ^(w) 

fi{z,z)a{w,w) -= — (B3d) 

y2\z — w] 1 / 4 

4>{z)o-{w,w)~ — — (B3e) 

y2(z — w) l ' z 

ifj{z)fi{w,w) ~ -=- — (B3f) 

y2(z — w) l i z 

J(-\ ( -\ Pl*P( W >™) /DO \ 

i/>(z)a{w,w) ~ —— — - (B3g) 

y2(z — w) 1 ' z 

T(-\ ( -\ Pl<r(w,w) , 

i/>{z)»{w,w) ~ —- — - (B3h) 

\J2[z — w) l i z 
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where some arbitrariness remains in the constants (3 and 7 because of the nonlocal character 
of these OPE's: (3 = ±1 and 7 = ±e ±Mr / 4 . In this work we choose 7 = e l7r//4 and (3=1. 

Since the regularized product of a with /i is a fermion, these operators must carry some 
anticommuting character. We shall assume that the disorder operator \i anticommutes with 
all fermion fields and other disorder operators, but not with itself nor with the order fields. 
This is a matter of convention [a could have been chosen instead). 
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FIGURES 
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FIG. 1. The 'zig-zag' chain, with interchain coupling k\ and intrachain coupling K2, also 

equivalent to a single chain with NN coupling k\ and NNN coupling K2. 
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